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Context. Elemental correlations and anti-correlations are known to be present in globular clusters 
(GCs) owing to pollution by CNO cycled gas. Because of its fragility Li is destroyed at the 
temperature at which the CNO cycling occurs, and this makes Li a crucial study for the nature of 
the contaminating stars. 

Aims. We observed 112 un-evolved stars at the Turnoff of the NGC 6752 cluster with FLAMES 
at the VLT to investigate the presence and the extent of a Li-0 correlation. This correlation is 
expected if there is a simple pollution scenario. 

Methods. Li (670.8 nm) and O triplet (771 nm) abundances are derived in NLTE. All stars belong 
to a very narrow region of the color-magnitude diagram, so they have similar stellar parameters 
(J eff , logg). 

Results. We find that O and Li correlate, with a high statistical significance that confirms the 
early results for this cluster. At first glance this is what is expected if a simple pollution of pristine 
gas with CNO cycled gas (O-poor , Li-poor ) occurred. The slope of the relationship, however, 
is about 0.4, and differs from unity by over 7<r. A slope of one is the value predicted for a pure 
contamination model. 

Conclusions. We confirm an extended Li-0 correlation in non evolved stars of NGC 6752. At 
the same time the characteristic of the correlation shows that a simple pollution scenario is not 
sufficient to explain the observations. Within this scenario the contaminant gas must have been 
enriched in Li. This would rule out massive stars as main polluters, and favor the hypothesis that 
the polluting gas was enriched by intermediate or high-mass AGB stars, unless the former can 
be shown to be able to produce Li. According to our observations, the fraction of polluting gas 
contained in the stars observed is a considerable fraction of the stellar mass of the cluster. 
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1. Introduction 

Among the light elements Li is special, because Li is produced in the Big Bang nucleosynthesis, 
enriched during the Galaxy evolution and destroyed in the stellar interior. 

After the pioneering works with 4m class telescopes (Pasquini and Molaro 1996, 1997, Castilho 
et al. 2000) , the study of Li in globular clusters (GC) has received a tremendous boost with 8m 



telescopes (IBoesgaard et al. 



1998 



Thevenin et al 



2001 



Bonifacio et al. 2002, 2007, Pasquini 
et al. 2004), and, with the advent of multi object s pectrographs, many un-evolved an d evolved 



stars have been observed for Li in the nearest GCs ( Korn et al. 2007 



Hernandez et al. 



2009; 



D'Orazi & Marino 2010; 



o2 



D'Orazi et al. 2010; 



Lind et al. 2009; Gonzalez 



Monaco et al. 2010). 



This work shall concentrate on the GC chemical inhomogeneities. Chemical inhomogeneities 
in globular clusters have been known for almost 50 years (Harding 1962) and Li may represent a 
key element in distinguishing between the possible mechanisms causing them. 

Historically, the first NTT observations of two Turn-Off (TO) stars in the metal rich, chemically 
inhomogeneous GC 47 Tuc did not show any difference in their Li content, even though one star was 
CN rich and the other CN poor (Pasquini and Molaro 1997). An early claim made by Boesgaard 
et al. (1998) of Li inhomogeneity in M92 has been questioned on the basis of a more realistic 
treatment of the errors (Bonifacio 2002). This situation has dramatically changed after the analysis 
of nine TO stars in NGC 6752, because for the first time strong Li variations were detected among 
these stars, and in addition they were anti-correlating with Na abundances; for seven stars an O 
measurement was also possible, showing a possible O-Li correlation (Pasquini et al. 2005) . 

It has been firmly established that chemical inhomogeneities in GC bring the signatures of 
CNO cycling, but all the possible reactions that involve the CNO cycle occur at temperatures well 
above 10 million degrees. At these temperatures, proton captures have destroyed all Li; the reaction 
destroying Li is indeed already very effective at 2.5 MK and has an extremely high power with 
temperature. The discovery of Li variations in NGC 6752 and their possible anticorrelation with 
processed material, but with the important detail that Li is still detected in the stars with the lowest 
O abundances, therefore suggests that the stars have been enriched by processed gas, possibly 
produced by the AGB stars in their H-burning shell phase. At these temperatures also the full CNO 
cycling takes place, and ox ygen is expected to be burnt. Indeed, an o xygen-sodium anticorrelation 



is observed in this cluster (iGratton et al. 



2001 



Carretta et al 



2007). We note that the NGC 6752 



TO stars have metallicity ([Fe/H]~-1.4) and effective t emperatures (Teff~ 6300 K) hot enough that 
they should be considered bona fide Spite plateau stars( Spite & Spite 1l982L for which internal Li 
depletion is excluded on the basis of standard models. If we were able to establish on a quantitative 
basis the existence of the expected Li-O correlation, we could quantitatively test the predictions 
of the main existing scenarios: these stars are the components of a second generation of low-mass 
stars in the cluster, which formed from the ejecta either of intermediate mass stars expelled during 
the AGB phase, as proposed by D' Antona et al. (2002), or by material processed by massive stars 
(Decressin et al. 2007). 

We are, on the other hand, conscious that not all clusters may have gone through the same 
formation process. NGC 6397 is showing so far a puzzling behavior: in this cluster Li is almost 
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constant in the observed TO stars and at the Spite plateau level, e xcept for a mild decrease with 
decreasing effective temperature ([Gonzalez Hernandez et alj|2009l) . while some stars show up to 
0.6 dex difference in their O abundance (Pasquini et al. 2004). O-poor-Li-rich stars have a puzzle in 
the pollution scenario, and only the production during the AGB phase would at least qualitatively 
explain it. 

In a simple self-pollution scheme Li is a powerful discriminator, and it becomes even more 
powerful when used in conjunction with O, (rather than Na or N) because both elements are de- 
stroyed in the full CNO cycle. For a given fraction of destroyed O at least the same fraction of Li 
must have been destroyed. 



2. Observations and analysis 



The observations were obtained with FLAMES and the GIRAFFE spectrograph located at the VLT 
Kueyen (Pasquini et al. 2003). The Li observations were obtained in ESO Period 73 and the O IR 
triplet observations in ESO period 75 for the same stars. The resolving power and S/N ratio obtained 
were of R ~ 17, 000, S/N ~ 80 - 100 for the Li (607.8 nm) region, and for O I (777.1-777.5 nm), 
R~ 18,400andS/JV~40-50. 

Targets were selected on the basis of their magnitudes and color and on their distance from 
the cluster center. T he BV photometry is based on the data collected in the framework of the EIS 
Pre-Flames survey (Mo manv et al] 1200 II) . The magnitudes were extracted with DAOPHOT and 
calibrated against the Stetson photometric sequence present in the same field. The following se- 
lection criteria were used photometrically: stars must be in the TO region with a magnitude range 
17.0 < V < 17.3 and color range 0.42 < ( B - V) < 0.46; spatially, in a radius 4 < r < 9 arcmin 
from the center of the cluster in order to limit the effect of crowding and background light con- 
tamination. Furthermore, each selected star has been visually inspected to guarantee the absence 
of nearby bright objects within 5 arcsec. With this restricted limit in magnitudes and colors the 
physical parameters (T e //, logg) of the targets are very similar. This shall ensure that the large 
chemical abundance differences observed cannot be attributed to different stellar parameters. The 
color-magnitude diagram of NGC 6752 highlighting the selected stars is shown in Fig. 1 . The stan- 
dard deviation in (B - V) color of our sample is 0.008 magnitu des, consiste nt with our error estimate 
of 0.01 mag. We adopted a reddening E(B — V) = 0.04 mag karris! 1 99^ and computed ef fective 



Gonzalez H ernandez & Bonifacio! (120091) . The 



temperatures from the IRFM based calibration of < 
mean effective temperature of our sample is 6353 K with a standard deviation of 38 K. As discussed 
above, this is entirely accounted for by the photometric error. We thus assumed all stars to have this 
temperature. For the surface gravity we assumed log g = 4.0. 

Li and O triplet equivalent widths were then measur ed by fitting Gaussia n profiles. The Li 

(2010), which is a result 



Sbordone et al. 



Ludwig et al. 



(2009) 



abundances were computed with the fitting function of 
of the NLTE synthesis and 3D hydrodynamical models of the CIFIST grid 
Unfortunately a similar fitting function for th e O triplet is not yet available. W e therefore computed 
the NLTE O abundances with the Kiel code ( Steenbock & Holwegerl 1984 ), a model atmosphere , 
computed with the Linux version ( Sbordone et al kool of the ATL AS 9 codeJKuruczl 19931 



no overshooting and the 1 km s 1 opacity dis t ribut ions functions of 



Castelli & Kurucz 



2005), 



the model atom described in 



Paunzen et al 



(2003), and 



(1999). The collisions with neutral H are treated in 
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an approximate way 
as recommended by 



(see 



Steenbock & Holweger 1984), and we adopted a scaling factor of 1/3 



Caffau et al 



(2008). The only Oi line that could be safely measured was the 
777.1 nm line, the strongest of the triplet. If the other lines of the triplet could also be measured 
they implied oxygen abundances that could disagree by as much as 0.5 dex from the strongest line. 
We therefore believe it is safer and mor e homogeneous to d erive oxygen abundances only from this 



line. We adopt the solar abundances of 



Caffau et al 



(2010). 



It has long been known that the behavior of the O i permitted triplet lines is a challenge to 
our comprehension. For example, Spite & Spite (1991) compared the oxygen abundances derived 
from the triplet and from the forbidden [O i] lines in LTE, finding discrepancies between 0.35 dex 
and 0.75 dex. We cannot use the [O i] lines because they are not covered by our spectra. And the 
lines should have equivalent widths (EWs) of less than 0. 1 pm, which is beyond the measurement 
capabilities for stars of this magnitude with VLT. However, Garcia Perez et al. (2006) have studied 
in detail a sample of bright subgiant stars and compared the oxygen abundances derived from 
different indicators, including the Oi triplet, computed in NLTE, and the [Oi]. They found on 
average an offset of 0.19 dex, the triplet yielding a higher abundance. It is likely that other effects, 
such as granulation, are the cause of this discrepancy, and work is in progress to address this with 
C05BOLD hydrodynamical models (Ludwig et al. 2009). For the purpose of this letter we stress 
that only the relative oxygen abundances are relevant. The narrow range in effective temperature 
and surface gravity covered by our sample strongly argues in favor of the notion that the observed 
variations in EWs of the O I lines are caused by oxygen abundance rather than to other effects. 

The errors on both Li and O are dominated by errors in effective temperature and in equivalent 
width measurement, because the errors due to surface gravity and microturbulence are negligible 
with respect to the former. We adopt an error of 100 K on the effective temperatures, which corre- 
sponds roughly to the semi-dispersion in the effective temperature s determined for our sample. The 
errors on equivalent widths were estimated with Cayrel's formula ( Cavrelfl988 ). The total error on 
abundance is obtained by summing under quadrature the errors associated to effective temperature 
and equivalent width. These amount to 0. 14 dex for oxygen and 0.09 dex for lithium (with opposite 
sign for the two elements). 

Our results are displayed in Fig.|2j the stars for which both oxygen and lithium could be mea- 
sured are shown as filled symbols. Open symbols with arrows are for the stars for which one or 
both measures are upper limits. We have 77 measurements of both elements, for 25 stars lithium is 
measured, but oxygen is not, for 4 stars oxygen is measured and lithium is not, finally for 6 stars 
neither element could be measured. This makes a total of 1 12 measurements and upper limits. 

The visual impr ession is that of a real correlation between the abundances of oxygen and 
lithium, as found by IPasquini et al. I (|2005), albeit in presence of a large scatter. This is indeed 
confirmed by the statistical analysis of the data. If we confine the sample to the 77 stars for which 
both elements are measured, the non-parametric test of Kendall's r tells us that there is a 98% 
probability that the two are positively correlated. To see if the upper limi ts may also bring som e 
information we used survival statistics, a s provided by asurv Rev 1.20 jLavallev et al. ( 1992 )). 
The generalized version of Kendall's t (Br own et al. 1119741) . as described by Isobe et al. 



(1986) 
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Fig. 1. Color-magnitude diagram of NGC 6752. The target stars are highlighted. Note the restricted 
range in magnitudes. 



provides a probability of correlation of 99.7%, reflecting the fact that the upper limits are mostly 
found along the correlation at low values of A(Li) and [O/Fe]. 

After establishing that there exists a correlation between [O/Fe] and A(Li), one may venture to 
perform parametric fits. A least-square fit using errors in both variables provides a slope of 0.40 
with an error of 0.08 and an intercept of 2.37 dex with an error of 0.01 dex. This slope corresponds 
roughly to what would have been guessed by eye and is shown as a dotted line in Fig. [2] What 
is interesting is that the root mean square deviation about this relation is 0.12 dex, which is even 
slightly lower than the sum under quadrature of the errors on oxygen and lithium. This strongly 
suggests that in reality there is a very tight relation between [O/Fe] and A(Li) and that the observed 
scatter is only due to observational errors. 

The slope of the relationship is important because, as mentioned in the introduction, the simple 
pollution scheme by CNO processed gas implies that Li and O shall diminish at least at the same 
rate. Therefore the slope should be one. This is clearly not the case, and the formal slope found by 
our fit is smaller, in fact a slope of unity is at over 7 <x from our best-fit value, therefore the case for 
a slope of one can be ruled out with very high confidence. We stress that this result does not depend 
on the unique effective temperature adopted for all stars; the same result is found when adopting a 
different temperature scale (based on Ha wings, Shen et al. 2010, in preparation). 



3. Discussion 



We can confirm in this study a Li-O correlation for NGC 6752. Such a neat correlation (togethe r 



with the Li-Na anticorrelati on) is not observed in NGC 6397 ([Fe/H]= -2.03 



Gratton et al. 



and47Tuc ([Fe/H]= -0.64; 



Carretta et al 



(2001)) 



(2004)), the other two clusters observed in a similar way. 



The correlation of NGC 6752 is close to what one would expect if a pollution scheme was at work. 



5 



2.6 - 



2.4 - 



d. 2.2 
< 



2.0 



1.8 - 



Z. Shen et al.: Li and O in NGC 6752 

i — i — i — i — i — | — i — i — i — i — | — i — i — i — i — r 




3 



<-Q 



"1 



J i i i i I i i i i I i i i i L 



-1 



-0.5 
[O/Fe] 



0.5 



Fig. 2. Li vs. O abundance for the observed stars. 



If O decreases by 0.8 dex without Li - enriched polluting material, then Li must decrease by at 
least the same quantity. Because we observe less Li depletion than O depletion, some Li must be 
produced by the contaminating stars, i.e. the polluting material cannot be Li-free. 

Another evidence for this con clusion is given by the observations of O-poor - Li-rich s tars in 



NGC 6397 dPasquini et al. 



2004). The most recent FLAMES results dD'Orazi et al 



2010) show 



very large variations of Na and O for stars with the same Li abundances, pointing again to a pollut- 
ing material with Li content. Our sample also shows in Fig.[2]a fairly large spread of O abundance 
for the same Li and suggests the presence of O-poor - Li-rich stars, but only for one star the devia- 
tion from the linear relationship is statistically significant. If the presence of O-poor - Li-rich stars 
were to be confirmed by higher S/N ratio observations, this would provide strong constraints on the 
maximum Li abundance in the ejecta. It is interesting to note that the Li enrichment is in general 
not so extreme as to destroy the general relationship, but the polluting gas should not have been 
homogeneous and well mixed in the whole cluster. Not all stars are in fact affected in the same 
way: only for a considerable fraction of them (5-10% in our sample), on the other hand, the gas 
polluting these stars should have been very Li-rich. If we consider a star depleted by 0.7 dex in O, 
this implies a percentage of polluted gas of 80% and only 20% of pristine one. 

The star with the highest lithium abundance has a value that is consistent with the prediction of 
the standard Big Bang nucleosynthesis and the measurement of the baryonic density of the WMAP 
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satellite (Cvburt et al. 2008). We believe this to be fortuitous, because we are invoking Li-polluted 
material to have contributed to the atmospheres of the studied stars and we are also stating that 
nothing can be firmly said on the "pristine" abundances of the elements affected by pollution. 
Although in general we may expect the Li content to be lower than the original value, there is no 
strong reason to say that, occasionally, it may not be higher, at least until the properties of the 
polluting stars are firmly established. 

One problem pointed out in several occasions for the formation from AGB ejecta is that the 
amount of 'polluted' material required is huge, and from our sample we can confirm that from 
the oxygen data of NGC 6752 the amount of polluted gas is comparable to half of the gas that 
formed the observed stars (Shen et al. 2010, in preparation). This can be seen in Figure 2: in this 
scheme, all stars with [O/Fe]<~0.0 are composed by at least half of the ejected material. This would 
pose serious constraints on the cluster IMF. Our sample should be unbiased as far as pollution is 
concerned, so it should represent the average chemical distribution of the cluster. 

We believe that any single pollution-dilution scheme we assumed could be too simplistic, and 
that the formation scenario was likely more complex. 

Acknowledgements. Z.X.S. acknowledges Chinese NSFC Grant 10903012. 
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